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A SIMPLE PREDICTION ON THE SURFACE DENSITY OF GALAXIES AT Z ~ 6 
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ABSTRACT 

Systematic surveys are being proposed to discover a significant number of gaiaxies at 2: ~ 6, which 
is now suggested as the epoch when the reionization era of the Universe ends. To plan such surveys, 
we need a reasonable expectation of the surface density of high redshift galaxies at different flux limits. 
Here we present a simple prediction of the surface density of 5.5 < z < 6.5 galaxies in the optical regime, 
extrapolating from what is already known about galaxies at 2 ~ 3. This prediction is consistent with 
the results of nearly all known searches for objects at z ~ 6, giving confidence that we may use it to plan 
optimal combination of survey depth and sky coverage in searching for such objects. We suggest that the 
most efficient strategy with existing ground-based facilities is to do medium-depth (m^s — 24.0 - 24.5 
mag) , wide- field (a couple of square degrees) survey using a wide-field camera at a 4m-class telescope. 
As the predicted surface density at this brightness level is very sensitive to the value of L*, the result of 
such a survey can be easily used to constrain the luminosity evolution from z ~ 3 to 6. 

Subject headings: galaxies: high-redshift — galaxies: luminosity function 



o 



X 



1. INTRODUCTION 

In the last several years, our knowledge about the Uni- 
verse at high redshift has been gradually extended to 
z ~ 6. As of today, five galaxies at z > 5.5 (Weymann 
etal. 1998; Hu etal. 1999, 2002; Dawson etal. 2001) and 
four quasars at z > 5.5 (Fan etal. 2000, 2001) have been 
spectroscopically confirmed. The complete Gunn-Peterson 
trough detected in the z = 6.28 SDSS quasar by Becker 
etal. (2001) and further investigation (Fan et al. 2002) led 
these authors to tentatively identify the end of reioniza- 
tion epoch of the Universe at z ~ 6. Thus the assessment 
of z ~ 6 galaxy number counts at different brightness lev- 
els will have a very direct cosmological impact, since it 
will quantify the number density of UV-emitting objects 
that are the physical cause of the reionization. Several sys- 
tematic surveys aimed at discovering a significant number 
of z ~ 6 galaxies have now been proposed and are be- 
ing carried out. To make such surveys efficient, we need 
to have a rough idea of the surface density of galaxies at 
z ~ 6. There are a few theoretical predictions on the sur- 
face density of galaxies at high redshifts, either based on 
N-body simulations (cf. Weinberg etal. 1999, 2002) or 
based on semi-analytic formalisms (cf. Robinson & Silk 
2000). However, those predictions are more qualitative 
than quantitative at this point, because they either are 
limited by finite volume and finite resolution, or have to 
rely on several parameters which remain very uncertain in 
the absence of significant amounts of actual data. Thus 



one may prefer not to base survey plans directly on such 
predictions at the moment. 

In this paper we present a simple, observational ap- 
proach. We may assume a reasonable luminosity function 
for galaxies at z ~ 6 by extrapolating the known results 
at z ~ 3, which is the highest redshift where the lumi- 
nosity function of galaxies has been quantified over a wide 
enough brightness range (Steidel etal. 1999). We will use 
existing data at z ^ 5 to constrain the normalization of this 
extrapolated luminosity function. Once the z ~ 6 lumi- 
nosity function is estimated in this way, the surface density 
can be calculated in a straightforward manner, i.e., by nu- 
merically integrating the luminosity function over the vol- 
ume occupied by unit sky-coverage in the redshift bin of 
interest. In §2, we present the details of such a surface 
density prediction over the range 5.5 < z < 6.5. We com- 
pare our prediction to all the available observations in §3. 
Comparison to two theoretical models is made in §4. A 
summary is given in §5. 

2. THE PREDICTED SURFACE DENSITY OF GALAXIES AT 

5.5 < z < 6.5 

2.1. The Extrapolated Luminosity Function 

We start with the z ~ 3 luminosity function from Stei- 
del etal. (1999), who merged the ground-based and the 
HDF-N Lyman-break galaxy samples. Since the z ~ 3 
sample used in Steidel etal. (1999) was selected based on 
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the Lyman-break signature in the SED of galaxies, our 
approach makes the imphcit assumption that galaxies at 
higher redshifts also manifest themselves by the Lyman- 
break feature. Because the line-of-siglit intervening H I 
absorption is more severe at higher redshifts and the ef- 
fective break moves from the 9 12 A Lyman limit to the 
Lya line at 1216A (cf. Madau 1995), one can expect that 
the Lyman-break signature is even stronger, therefore our 
assumption is likely valid. However, there might exist a 
different type of galaxy whose UV-pliotons are completely 
absorbed by dust, like the ultraluminous infrared galaxies 
that are known to exist at lower redshifts. There is no 
Lyman-break for such galaxies, since they have essentially 
zero flux from the UV to the optical range. Our predic- 
tions neglect these objects, becaiisc they are also absent 
from the Lyman-break samples at 2 ~ 3. Since we deal 
with the emitted UV, such objects will not affect our con- 
clusions. 

The details of the luminosity hmction extrapolation 
from 2; c:i 3 to z ~ 6 depend on the adopted cosmologi- 
cal model. Different cosmological parameters will change 
the apparent magnitude of an object with a given intrinsic 
brightness, and will also affect the physical volume cor- 
responding to a given sky coverage. We consider three 
cosmological models, namely, a flat universe without cos- 
mological constant (Qm = 1 and f^A = 0), a flat, low 
mass universe with J7m = 0.3 and J^a = 0.7, and an 
open universe without cosmological constant = 0.2 

and = 0). For simplicity, we denote the three models 
by their Om and Qa values as the (1, 0), (0.3, 0.7), and 
(0.2, 0) models, respectively. Throughout the paper we 
use a Hubble constant of Hq = 65 km s^^ Mpc~^. All the 
quoted magnitudes are AB magnitudes and, unless noted 
otherwise, refer to the 1400 - 1500A spectral range in the 
rest-frame. 

The luminosity function of z ~ 3 galaxies of Steidel 
etal. (1999; see their section 4.3) is a Schechter function 
of the form 4>(L) = (4>*/L*)(L/L*)"ea;p(-L/L*). The L* 
galaxies at this redshift have 7^-band apparent magnitude 
of m* = 24.48 mag. At z = 3.04 (the median redshift 
for Steidel et al.'s sample), this m* value corresponds to 
M* = -20.37, -21.23 and -21.31 mag in the (1, 0), (0.3, 
0.7) and (0.2, 0) models, respectively. The 7^-band has 
central wavelength of Aq = 6930A, or ^ 1700A in the rest- 
frame at z ~ 3. Since the SED of a young non-dusty 
galaxy in fi, is essentially flat from 1400 - 1700A, these 
magnitudes apply to the 1400 - 1500A range as well. As 
a first approximation, we also assume that there is no sig- 
nificant luminosity evolution for galaxies from z ~ 3 to 
5.5 < z < 6.5, such that L* galaxies at 5.5 < z < 6.5 still 
have the above mentioned absolute magnitudes. Possible 
evolution effects will be discussed in §5. We noticed that 
Steidel etal. (1999) made a similar assumption in com- 
paring the ground-based + the HDF-N sample at z ~ 4 
to its counterpart at z ~ 3, concluding that the observed 
z ~ 4 galaxy distribution was completely consistent with 
such an assumption at the bright end, and within a factor 
of two at the faint-end. In fact, they suggested that the 
faint end mismatch was due to the genuine structure in the 
small HDF-N field, rather than a true luminosity function 
discrepancy. We use continuum magnitudes, ignoring the 
possible contribution of Lya line emission, for two rea- 



sons: 1) only 50% of the galaxies in the sample of Steidel 
et al. are Lya emitters, although there was a strong bias 
in favor of such objects in the spectroscopic identification, 
and 2) the vast majority of those Lya emitters are weak 
in Lya strength. We assume a slope of a = —1.6, but our 
conclusions are insensitive to the exact value. 

Next we choose a normalization to fix the luminosity 
function, which amounts to allowing density evolution. 
Observationally, this can be done by adjusting $* such 
that the calculated surface density, either differential or 
cumulative, matches the availible observations. This is dif- 
ficult in our case, since there is no precise observed value 
to be used. There are only four published z > 5.5 galax- 
ies with spectroscopic information. These are the z=5.60 
galaxy in the HDF-N (Weymann etal. 1998), the z=5.74 
galaxy in the Hawaii Survey Field SSA22 (Hu et al. 1999), 
and the z=5.767 and z=5.631 galaxies in the HDF-N fiank- 
ing fields (Dawson etal. 2001). Statistically, these results 
are inadequate for the purpose of normalization, because 
their selection is not readily quantifiable. Nevertheless, we 
can still make a rough estimate of the cumulative surface 
densities of z > 5.5 galaxies in the well-studied HDF-N and 
use this value as our normalization. We consider two cases: 
a cumulative surface density of 1.37 per arcmin^ and 0.11 
per arcmin^, both to a limit tuab = 27.0 mag. The for- 
mer case (high normalization) is equivalent to one z > 5.5 
galaxy per NIC-3 field (cf. Thompson etal. 1999) and 
the later one (low normalization) is equivalent to one such 
galaxy per WFPC2 HDF coverage. Note that toab = 27.0 
mag is not the magnitude of this z=5.60 galaxy, but an es- 
timate of the selection limit for a galaxy at z > 5.5 as ob- 
served at 9100 - 9800A. This galaxy is not a Lya emitter, 
confirming that the magnitude in our prediction should be 
taken as referring to the continuum brightness level rather 
than including this emission line. 

2.2. The Predicted Surface Densities 

Cumulative surface densities can now be calculated by 
numerically integrating the luminosity function over the 
volume occupied by unit sky coverage in the range 5.5 < 
z < 6.5. Fig. 1 shows these results, where the cumulative 
surface density is presented as the total number of galaxies 
per deg^, whose brightness is brighter than a specified ap- 
parent AB magnitude. The predictions for the (1, 0), (0.3, 
0.7) and (0.2, 0) models are plotted in solid, long-dashed 
and short-dashed lines, respectively. Thick lines are used 
for the high normalization case, while thin lines are used 
for the low normalization one. The absolute magnitude 
scales arc labeled on top of the figure, from bottom to the 
top for the (1, 0), (0.3, 0.7) and (0.2, 0) models, respec- 
tively. The M* values for these models are labeled in the 
legend in the parenthesis together with the correspond- 
ing apparent magnitudes, and are also indicated along the 
absolute magnitude scales at the top by upward arrows. 
Numerical results are presented in Table 1, where the num- 
ber density is listed as number of galaxies per arcmin^ for 
straightforward use in planning future observations. 

3. CONSISTENCY CHECK FROM LIMITED OBSERVATIONS 

As a minimal consistency check, we compare our pre- 
dictions to the limited (direct and indirect) observations 
in hand for galaxies at z ~ 6. 
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Fig. 1. — The cumulative surface density predictions of 2 ~ 6 galaxies in the 5.5 < z < 6.5 redshift bin, presented as the total number of 
galaxies per deg'^ whose brightness is brighter than a specified apparent magnitude (defined around 9100 - 9800A on the AB scale). The 
results for the (1, 0), (0.3, 0.7) and (0.2, 0) models are plotted in solid lines, long-dashed lines and short-dashed lines, respectively. The 
high and low normalization cases are plotted in thick and thin lines, respectively. The "-h" labels mark the positions where the HDF-based 
normalizations are placed. The corresponding absolute magnitude scales are labeled on top of the figure, from bottom to the top for the (1, 
0), (0.3, 0.7) and (0.2, 0) models. The M* values in these three models are given in the legend together with the corresponding apparent 
magnitude scales, and are also marked with the upward arrows on the scales. The upper limit derived from the preliminary data of our survey 
is indicated by the thick downward arrow, while the lower limit derived from Hu et al. (1999) is labeled by the thick upward arrow. The 
indirect bright-end upper-limit inferred from the hosts of the SDSS 5.5 < z < 6.5 quasars is indicated by "±" from the left to right for (1, 0), 
(0.3, 0.7) and (0.2, 0) models, respectively. 



3.1. The Intermediate Brightness Level (23 < niAB < 24 

mag) 

Our group has started a multi-color search for 2 ~ 6 ob- 
ject, aimed at constraining the luminosity function at this 
intermediate brightness level. Four medium-band filters, 
namely, m(802nm), n(848nTO), o(919nm) and ]3(974nm), 
are used in this survey. These are the four reddest of 
our 15 medium-band filters which cover the entire 300 - 
1000 nm range, and arc designed to avoid the brightest 
and most variable night-sky lines to yield dark background 
and fringe- free imaging. These four filters are ideal for 
discriminating the Lyman-break signature in the SED of 
galaxies at 5.5 < z < 6.5. Specifications of these filters 
can be found in the scries of papers of the BATC (Beijing- 
Arizona- Taiwan- Connecticut) Survey, whose filters arc es- 
sentially the same but smaller in physical size (Fan ct al. 
1996; Shang etal. 1998; Zheng etal. 1999; Yan etal. 
2000). 

Our medium-band data so far allow us to put a mean- 
ingful constraint on the luminosity function at z ^ 6. In 
our pilot run in June 2000, we observed a 35' x 35' area 
around the HDF-South in the reddest three bands (n, o 
and p), using the MOSAIC-II array mounted at the CTIO- 
4m telescope. Details of the observation will be presented 
elsewhere (H. Yan et al., in preparation). To summarize, 
we reached AB magnitudes of n(848nm) = 23.2, o(919nm) 
= 23.5 and p(974nm) = 23.0 mag at 5-cr with > 50% com- 
pleteness, and used these images to search for drop-outs at 
^; ~ 6. We found three candidates which were absent from 



the n(848nrn)-band image but prominent in the o(919nm)- 
band. However, since these objects are also very faint in 
the p(974nTO)-band, imaging in bluer passbands is needed 
to confirm their nature. Two out of these three objects 
lie within the broad band UBVRI CTI0-4m BTC images 
of Palunas etal. (2000), whose survey field has about a 
15' X 15' overlapping region with ours. These two drop- 
outs are both clearly visible on the BTC _R-band and short- 
ward images, but not on the I-band image. Thus these two 
objects are most likely lower redshift interlopers, possibly 
z = 2.3 quasars with Mgll redshifted to ~ 9100A. This 
meant that at best only one positive candidate remained 
that could be at z ~ 5.5 6.5. 

Based on these results, we can put a safe upper limit 
to the cumulative number density of 5.5 < z < 6.5 ob- 
jects as one object in 35' x 35' down to the flux limit of 
wiAB = 23.5 mag. This upper limit is indicated on Fig. 1 
by an downward arrow, and is consistent with both high 
and low normalization cases in our prediction. 

3.2. The Faint End {niAB > 25 mag) 

A lower limit at this brightness level can be obtained 
from the z = 5.74 galaxy of Hu etal. (1999), assuming 
that its single strong emission line identification is reliable. 
The galaxy's Z-band magnitude (free of the emission line) , 
which is similar to the SDSS z'-band magnitude, is 25.5 
mag. Since their survey area is about 390" x 280", a simple 
lower limit estimation is log{N) = 2.07/ deg^ to m = 25.5 
mag. This lower limit is indicated on Fig. 1 as "_L", and 
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fits with our high normalization case but not with the low 
one. 

Unfortunately, the two z > 5.5 galaxies of Dawson et al. 
(2001) cannot be used as a direct constraint to the lumi- 
nosity function, since no continuum magnitude was given 
for either of the two. However, because their Iab magni- 
tudes are all fainter than 25, their result is at least not in 
conflict with our estimates above. 

3.3. Indirect Result at the Bright End {ruAB ~ 23 mag) 

Finding any z ~& galaxies at these bright levels would 
require a SDSS-like all-sky survey, but extending at least 
one magnitude deeper, which is a daunting, if not impos- 
sible, task with any existing facilities. However, we can 
get some information in this regime by using quasar host 
galaxies as the tracers of the entire galaxy population. To 
the extent that quasar hosts can be taken as represen- 
tatives of galaxies, it is possible to obtain a reasonable, 
although indirect, bright-end lower limit if the hosts lumi- 
nosity of the known z ~ 6 quasars can be obtained. It has 
been long suggested that there is a linear relation between 
the quasar luminosity and the minimum luminosity of its 
host galaxy (cf. McLeod & Rieke 1995), which can be un- 
derstood in the sense that a more luminous host galaxy 
is required to fuel a more luminous quasar. Thus there 
might be a similar relation between the luminosity of the 
quasar and that of its host at higher rcdshift. Such a rela- 
tion would give a luminosity function constraint from the 
four z > 5.5 quasars known, all discovered by SDSS in ~ 
1500 deg2 of survey area (Fan etal. 2000, 2001). 

We construct such a relation from Bahcall et al. (1996), 
who imaged twenty nearby (z < 0.3) luminous quasars 
with HST, and obtained their host luminosities. A sim- 
ple fit to their data gives Mhost — 0.425 x Mqso — 11-82 
{rms = 0.38 mag). Using this formula, we found that 
the host luminosities of the four z > 5.5 SDSS quasars 
were almost identical. Their average absolute magnitude 
at around rest-frame 1400A is M = —23.1, —23.5 and 
-23.7 mag (rms ~ 0.2 mag) in the (1, 0), (0.3, 0.7) and 
(0.2, 0) universe, respectively, which correspond to appar- 
ent magnitudes of 22.8, 23.4 and 23.6 mag, respectively, 
at z = 6. The AGN unification scheme implies that if an 
AGN is detected, statistically there must be two or more 
galaxies which also harbor AGN activity, which cannot 
be detected because obscuring matter in the surround- 
ing torus largely blocks our view. In other words, since 
there are four z > 5.5 quasars detected in ^ 1500 deg^, 
the total number of z > 5.5 galaxies (including the four 
quasar hosts) is at least 12, implying a number density 
of log{N) > 2.10/ deg^ at the above mentioned brightness 
levels. As indicated above, this number should be taken 
as a lower limit, because there could be some galaxies as 
luminous as the quasar hosts, but lack a massive central 
black hole and thus cannot be traced by AGN. This in- 
direct constraint is shown by the upward arrow in Figure 
1. 

There are two major uncertainties on this bright-end 

constraint. One is the effect of quasar luminosity evolu- 
tion. The above linear relation between the hosts and the 
quasars is inferred from a local sample, yet we apply it 
to z ~ 6. Typical quasar luminosities have been shown 
to evolve with redshift, as (1 -|- z)^ with /3 = 2.5-4 for 



z < 2.5. This indicates that quasars may have once been 
brighter with respect to their hosts than we see today. In 
other words, the hosts of the z ~ 6 quasars could be fainter 
than we estimated here, and therefore the limits shown in 
Fig. 1 could be moved more to the right. However, since 
it is almost impossible to obtain a quantitative estimate 
about the amplitude of such a shift, we just leave the de- 
rived constraint as it is now. Another issue is whether 
quasar hosts are fair representatives of field galaxies, and 
therefore whether the above bright-end limit is meaning- 
ful. Bahcall et al. (1996) speculated that the quasar hosts 
might not fit a Schechter function and might be 2.2 mag 
brighter than field galaxies. However, their suggestion was 
based on a volume-limited sample, which is clearly not 
the case for their luminous quasar sample. Furthermore, 
galaxy evolution clearly cannot be neglected in comparing 
results from 2;«0.2to2:~6. Net evolution of the hosts 
would move this constraint to higher luminosity, running 
counter to the effect of quasar luminosity evolution. 

Even in the light of these caveats, we believe that the 
bright-end constraint indicated by the SDSS z > 5.5 
quasars is useful. In fact, the (1,0) model, either high or 
low normalization, and the low normalization case of the 
(0.2, 0) model, are clearly not consistent with this con- 
straint. The high normalization case of both the (0.3, 0.7) 
and (0.2, 0) models, on the other hand, are perfectly con- 
sistent with this constraint. The low normalization case 
of the (0.3, 0.7) model is barely consistent with this con- 
straint. Thus, the mere detection of four quasars at z > 5.5 
implies a rather high normalization to the counts and lu- 
minosity function of galaxies at these redshifts. 

3.4. Recent Narrow-band Search 

The only reported data which seem hard to reconcile 
with our prediction are from the z = 5.7 Lya emit- 
ter search of Rhoads & Malhotra (2001), who used two 
narrow-band filters to look for Lya emission. They re- 
ported 18 robust 2;=5.7 Lya emitters in 0.36 deg^ to a 5-0- 
survey limit of mAB=24.8 mag (average for both bands), 
and discussed cosmological implications of this result. Al- 
though these objects have not yet been spectroscopically 
confirmed, we compared our prediction against this result 
for the sake of completeness. For simplicity, we only com- 
pare it against our high normalization (0.3, 0.7) model. 

To make such a comparison, we need both their cumu- 
lative number density corrected for the survey volume and 
the 9100-9800A continuum magnitudes of their objects. 
The two narrow-band filters that they used, 815nTO and 
823nTO, are both 75A wide and can only pick up Lya emis- 
sion from z = 5.70 5.78. The co- moving volume of their 
survey is 3.27 x 10^ Mpc^ for their 0.36 deg^ sky coverage. 
Our model uses 5.5 < z < 6.5, giving a co-moving volume 
of 3.96 X 10^ Mpc^ for the same sky coverage. Since their 
survey volume is only 8% of this, the 18 objects that they 
detected should only account for 8% of the total galaxies in 
our prediction. In other words, if the 18 objects were taken 
at their face values, the total number of galaxies that our 
model would predict is 225 in our volume at 5.5 < z < 6.5. 
Furthermore, the objects they found are believed to be 
strong Lya emitters. Although spectroscopic identifica- 
tion is strongly biased in favor of such objects, past studies 
indicate that they are only a minority of the entire popu- 
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lation. If the 2: ~ 3 Lyman-break galaxy result of Steidel 
etal. (1999) is universal only ^ 25% of the Lyman-break 
galaxies have Lya emission strong enough to be picked up 
in such narrow-band searches. This means the total num- 
ber of z ~ 5.7 galaxies inferred from Rhoads & Malhotra's 
result is about 72, and thus the total number of objects 
that our model would predict is 900. Translated into num- 
ber density, it amoimts to 0.694 per arcmin^. Note that 
we assume their 18 objects form a complete sample down 
to their detection limit, which is possibly not the case. If 
we made a correction for their survey incompleteness, the 
inferred number density would be even higher. 

The continuum magnitudes of these objects were not 
explicitly given in their paper, therefore we use the equiv- 
alent width of the weakest object to estimate their survey 
limit in terms of continuum magnitude. Using the numbers 
from Figure 1 of Rhoads & Malhotra (2001), we find that 
the weakest emitter has fi,{line) = 0.48 /iJiy in the line, 
or equivalently mAB{line) = 24.70 mag, almost approach- 
ing their detection limit. Using the approximation that 
fv {line) / f I, {continuum) = {W/AX) + 1, where W is the 
equivalent width and AA is the bandwidth, we find that 
f I, {continuum) = 0.23/iJjy, or m^B {continuum) = 25.50 
mag. At this brightness level, our model gives the cumu- 
lative surface density of only 0.138 per arcmin^. 

This comparison shows that Rhoads & Malhotra's re- 
sult gives a factor of five times higher cumulative number 
density as our prediction does. Stated differently, if their 
result were used as our normalization at m^s = 25.50 
mag, the number of z 6 objects in the HDF-N would 
be much higher (> 6 times) than what have been found. 
Since the Rhoads & Malhotra's objects still need spectro- 
scopic identification to judge their real nature, we shall not 
go further than pointing out that there is a discrepancy. 

4. COMPARED WITH THEORETICAL MODELS 

We compare our estimates with the hierarchical N-body 
modeling results of Weinberg et al. (2002) in Fig. 2. The 
thick solid line is our (0.3, 0.7) high normalization esti- 
mate. The thin lines of different types are the predictions 
of Weinberg etal. (2002), reproduced by reading off the 
numbers from their Figure 8 and following the convert- 
ing procedures given by their paper. Only three of their 
models are reproduced here, namely, CCDM, LCDM and 
OCDM, which correspond to (f7M,f^A) = (1,0), (0.4,0.6) 
and (0.4,0), respectively. Their choices of 51m and Ha 
are slightly different from ours, but the effects of these 
on the global trends are only marginal. Two flavors of 
these models are shown: one is without any dust extinc- 
tion {A = at restframe 1500 A), and the other is with 20% 
dust extinction {A — 0.2 mag). Note that the count pre- 
dictions can differ widely solely due to different assump- 
tions about the internal extinction. Obviously only the 
(0.4, 0.6) model (LCDM) does not conflict with the HDF- 
N result as summrized in our Fig. 1 and §2.1. Further- 
more, nearly all those models predict unrealistically high 
surface density at mAB< 24.5 mag, which is caused by 
the finite resolution in the simulations. Trading volume 
for finer resolution (Weinberg etal. 1999), their LCDM 
result would give a more realistic surface density at the 
bright end, but again would predict too high a value at 
the faint end. 



It is also interesting to do a comparison among the the- 
oretical predictions. For example, we can compare the 
result derived from Weinberg etal. (2002) as shown in 
Figure 2 with the semi-analytic result of Robinson & Silk 
(2000) as shown in their Figure 3. A direct, quantita- 
tive comparison is not meaningful, since the result here is 
for 2; = 6 while the highest redshift given in Robinson & 
Silk's Figure 3 is z = 5, and their choices of Qm and Qa 
are slightly different. But we can still compare the global 
trends, which show obvious differences. In Robinson & 
Silk's semi-analytic formalism, the (1, 0) model gives the 
lowest galaxy counts, the (0.3, 0.7) lies in between, and 
the (0.2, 0) one gives the highest galaxy counts. They 
explained such trends as the result caused by the domi- 
nant effects of the different volume element dV/ dz and the 
growth factor G in different cosmologies. A lower-density 
universe always has higher dV/dz and larger G, and the 
combination of these two factors dominates the compet- 
ing effect of the larger luminosity distance which makes 
objects fainter in such a universe. On the other hand, 
Weinberg et al.'s (1, 0) model yields the highest counts, 
(0.4, 0.6) yields the lowest counts, and (0.4, 0) is in be- 
tween. Notice that in our approach (1, 0) gives the high- 
est counts, (0.3, 0.7) is in between, and (0.2, 0) gives the 
lowest counts. It is possible that the reasoning of Robin- 
son & Silk (2000) might only be partially true, and that 
other effects which they did not consider might actually 
be important. For example;, if the luminosity of L* galax- 
ies in a high-density universe is brighter than that can be 
produced from their current model, the number counts in 
different cosmologies might have a different behavior from 
that shown in their paper. 

5. DISCUSSION AND CONCLUSION 

5.1. The Effect of Luminosity Evolution 

The most crucial assumption made in our prediction is 
that the luminosity evolution from z ~ 3 to 6 is not signif- 
icant, so that the value of L* at around restframe 1400A 
is still the same at z ~ 6 as at ~ 3. However, there are 
several possibilities where this condition could break down. 
There are at least two major competing effects which are 
relevant. One possibility is that the merger/star- forming 
rate could be lower at 2; ~ 6, and so would bring down 
the value of L*. On the other hand, both dust extinction 
and metallicities could be lower at earlier epochs as well, 
which would make L* brighter. Since these effects tend to 
cancel each other out, we believe that to first order our as- 
sumption is reasonable. Needless to say, the reality could 
be more complicated that what we assumed. For exam- 
ple, it has been suggested that L* evolves mildly in the 
form L* (X (1 + z)^ , where /3 varies from —1 to —1.5 (e.g. 
Lanzetta etal. 1999). As described below, the surface 
density at the bright end is extremely sensitive to the L* 
value. Therefore, a direct comparison between this simple 
prediction, which serves as a first approximation, and any 
future observations can give quantitative estimates on how 
the luminosity evolves from 2 c± 3 to 6. For the sake of 
completeness, here we discuss how the possibilities men- 
tioned above would affect our prediction. 

The time interval between z = 3 and z = 6 is about 
1.28 Gyr in a Qm = 0.3, Qa = 0.7, and Hq = 65 universe, 
and is likely only sufficient to allow one major merger (c/. 
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Fig. 2. Comparing our estimates with the N-body simulation prediction of Weinberg et al. (2002). The high normalization case of our 
(0.3, 0.7) model is plotted as the thick solid line. The thin lines of different types are the predictions of Weinberg et al. (2002), reproduced 
by reading off the numbers from their Figure 8 and following the converting procedures given by their paper. Two flavors of these models are 
shown: one is without any dust extinction (A = at around f500 A), and the other has a dust extinction oi A = 0.2 mag. Note how large 
the differences in the number count predictions are with only slight differences in the assumed dust extinction alone. Obviously only the (0.4, 
0.6) model (LCDM) does not conflict with the HDF-N result. 



Makino & Hut 1997). Hence a higher merger rate alone 
would at most make L* at 2 ~ 3 twice as bright as at z ~ 6. 
On the other hand, a higher merger rate would certainly 
make star formation rate higher, and this would further 
contribute to the luminosity. The later effect, however, 
has not yet been well quantified. As a very rough esti- 
mation, the higher star formation rate could contribute 
another factor of few in increasing L* at 2; ~ 3. Thus 
the overall effect of merger/star formation rate difference 
would make L* at z ~ 3 four to six times as bright as at 
z ~ 6, or a difference of 1.5-2 mag. 

The dust content and metallicities are likely to increase 
from 2; ~ 6 to 3, and these would affect L* in the oppo- 
site way. For example, one of the more reddened object 
in Steidel et al.'s sample, MS 1512-cB58, is quoted hav- 
ing E{B -V) ~ 0.3 mag (Pettini etal. 2001). Assuming 
the extinction law of Calzetti etal. (2000), this number 
translates to A ° ~ 2.6 mag. This means the galaxies 
at z ~ 6 could get brighter by 2.6 mag at most, if dust 
extinction is largely absent at this redshift. The metal- 
licity of galaxies at large redshift is again very hard to 
quantify and seems to have a wide spread (c./. Nagamine 
etal. 2001); but in any case, it is not likely that this effect 
would contribute more than 0.5 mag increase in brightness 
at 2 ~ 6 at around rest-frame 1400 A. 

To investigate how the surface density could be affected 
by the luminosity evolution, we also calculated the surface 
densities for different M* values. Specifically, we looked 
at the limits of the M* value where the high normaliza- 
tion case starts to conflict with the known constraints, and 
whore the low normalization case begins to be consistent 
with those constraints. Since Qm = 0.3 and = 0.7 are 
currently the most widely accepted values, for the sake of 
simplicity, we will only discuss our (0.3, 0.7) model here. 
In the high normalization case, if M* is brighter by more 
than 0.7 mag at z ~ 6, the predicted counts will conflict 
with our CTIO upper limit. If M* is fainter by only 0.3 



mag at z ~ 6, on the other hand, the counts will conflict 
with the lower limit derived from the SDSS QSO hosts by 
more than a factor of two. In the low normalization case, 
M* needs to be brighter by at least 2.0 mag to make the 
counts consistent with the Keck lower limit. In the mean 
time, however, such M* value makes the counts inconsis- 
tent with our CTIO upper limit, overpredicting the counts 
by a factor of two. This situation further confirms that the 
low normalization case can be rejected. 

5.2. Summary 

We present a simple empirical approach to predict the 
galaxy surface density at z ~ 6, which extrapolates the 
known luminosity function of z ~ 3 galaxies to z ~ 6. 
Our approach is based on only two observational results, 
namely, the observed luminosity function of z ~ 3 Lyman- 
break galaxies and the number of 5.5 < z < 6.5 galaxies 
in the HDF-N, and the assumption that there is no strong 
luminosity evolution for galaxies from z ~ 3 to z ~ 6. The 
biggest uncertainty in our estimates comes from the nor- 
malization, i.e., the actual number density of z ~ 6 galax- 
ies in the HDF-N down to the limit of uiab = 27.0 mag, for 
which we used one per WFPC-2 field and one per NIC-3 
field as our low and high normalization, respectively. We 
checked our results against the constraints derived from 
all known observations. It seems that the low normaliza- 
tion case can be rejected. The high normalization case, on 
the other hand, is consistent with most constraints if the 
prediction is made in the (0.3, 0.7) or (1, 0) models. The 
only observation with which our predictions do not agree 
is the narrow-band z = 5.7 Lya emitter result reported 
by Rhoads & Malhotra (2001), whose number density is 
at least 5 times higher than our prediction. If their re- 
sult were used as the normalization, it would suggest that 
the number of z > 5.5 objects in the HDF be at least 
6 times as many as actually found. Since these emitters 
still need future spectroscopic identification to judge their 



real nature, we conclude that this is only a potential con- 
flict. On the other hand, a direct comparison between our 
prediction and any future observations can give quantita- 
tive estimation on how the luminosity evolves at diff'erent 
epochs, as indicated in the previous section. 

To summarize, we believe that our simple approach can 
be used to plan future surveys, where there will always be 
compromise between depth and sky coverage. As our pre- 
diction indicates, currently the most realistic way to find 
a significant number of z ~ 6 galaxies with the available 
ground-based facilities is to do multi-color, medium-depth 
and wide-field surveys reaching continuum itiab ~ 24.0 - 
24.5 mag from 8400A to the CCD Q.E. cut-off at around 
1/xm, and covering a couple of square degrees. There are 
now several wide-field CCD cameras available at telescopes 
of sufficient light-gathering power, e.g., the MOSAIC-I/II 
at the KPNO/CTIO 4m's, the CFH12K at the CFHT and 
the Suprime-Cam at the Subaru. Carefully designed sur- 
veys at a 4m class telescope could possibly discover a few 
dozen L > L* z ~ 6 galaxies within a few nights of ob- 
servation (Yan etal. 2002, in preparation). On the other 
hand, deep, pencil-beam surveys from the ground are not 
likely to be very successful even with 8-1 Dm class tele- 
scopes. As Table 1 indicated, pencil-beam surveys with 
a few square arcmin field of view would have to reach at 
least niAB = 27 mag in the difficult spectrum regime red- 
der than 8400A to discover a significant number of such 
objects. Since at least two bands of observation at sim- 
ilar depth are needed to select drop-out candidates, the 
telescope time required is very costly if not unrealistic. 
In the immediate future, it should be possible to use the 
the Advance Camera for Surveys, which was installed on 
board HST in Macrh 2002, for drop-out searches to better 
constrain the luminosity function at z ~ 6. Accessing the 
faint end of the luminosity function at 2; > 6 will be one 
of the major goals that we will pursue with the Next Gen- 
eration Space Telescope (NGST) out to redshifts as high 
2; ~ 9 — 10, and possibly laeyond. 
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Table 1 

Number Density Prediction for Galaxies at 5.5 <« < 6.5 
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Note. — Predictions are given for {Qm,^a) = (IjO); (0.3,0.7), and (0.2, 0) cosmologies, m is apparent magnitude and M is absolute 
magnitude in the corresponding cosmology model. S is cumulative number density in arcmin"^, calculated for the entire 5.5 < z < 6.5 bin. 



